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ABSTRACT 


The spectral reflectivity of clouds of ice spheres is 
calculated for the near infrared wavelength region ( I a £ A £ V A\ ) 
for comparison with laboratory observations of ice clouds. The 
laboratory cloud particles were highly nonspherical ice crystals, 
however, it is shown that the major feature of the theoretical 
spectrum are in close agreement with the observations when spheres 
of equal volume are employed in the radiative transfer problem. 

The results provide some support for the contention that ran- 
domly oriented nonspherical particles may be replaced by equiva- 
lent spheres for calculations of intensities. 

The reflectivities are shown to be sensitive to the 
particle size, especially the reflectivity minima near 1.5 and 
2.0 m • Changes by a factor of two in the particle size cause 
significant variations in the major reflectivity features? 
this verifies that spectral reflectivity measurements can be 
of considerable value for cloud identification. The agreement 
of the theoretical and experimental spectra also indicates that 
it should be possible to obtain meaningful conclusions from 
infrared reflectivity measurements of planetary atmospheres, such 
as the observations which have been made on the atmosphere of 


Venus . 
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INTRODUCTION 

In the wavelength region from 1 to 4 microns the absorbtivity 
of ice, as well as that of many other materials, has an absolute 
value such that the reflectivity of clouds composed of micron- 
sized particles should show large variations with wavelength 
and with particle size. Laboratory measurements confirm this and 
hence suggest that near-infrared reflectivities may be useful for 
the identification of clouds in planetary atmospheres; it is con- 
ceivable, for example, to derive the particle composition and size 
and the cloud optical thickness from the spectral reflectivity. 
Sagan and Pollack (1967) , and later Hansen and Cheyney (1968) , have 
attempted to extract such cloud parameters from a comparison of 
theoretical computations to atmospheric observations; however, to 
determine whether such a procedure can yield valid results it is 
important to find if the computations can in fact duplicate 
laboratory spectra in which the cloud parameters are known. 

The primary purpose of this paper is to show that the major 
features in the best existing laboratory spectra can be duplicated 
by theoretical calculations. The results also indicate the 
sensitivity of the near-infrared reflectivities to variations of 
the cloud parameters. 

Several authors have previously made computations for single 
scattering by cloud particles in the near-infrared. A two-stream 
approximation for the multiple scattering has been developed by 
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Sagan and Pollack (1967) and applied to the near-infrared 
reflectivity of Venus. Good observational spectra in the near- 
infrared have been published both for the atmosphere of Venus 
(Bottema et al., 1965; Kuiper, 1962) and for terrestrial clouds 
(Blau et al., 1966). 

LABORATORY OBSERVATIONS 

Zander (1966a, b) made laboratory measurements which are well 
suited for a comparison to theoretical computations. The experi- 
mental methods he employed and a detailed description of his 
results will be published elsewhere [Zander, 1968]; hence in this 
section only the experimental conditions which most directly 
affect the theoretical interpretation will be given. 

Zander measured the near-infrared spectral reflectivity 
of an optically dense ice cloud at an angle of incidence 
d 0 = 15° from the normal, an angle of reflection fl = 15°, and an 
azimuth angle difference co - cp 0 = 0°; i.e., the measurements were 
for the specular direction. Evidence of a specular component in 
the diffuse reflection from terrestrial clouds and hazes has been 
reported (Tank, 1964; Deirmendjian, 1968) and hence observations 
in that direction may be an undesirable special case. However in 
the laboratory measurements it is unlikely that there was 
significant reflection from flat surfaces of horizontally 



oriented crystals because the particles were small and 
continually tumbling about [Zander, 1966b]. 

Zander (1968 and personal communication) sampled the 
particle shape and the size distribution by inserting a 
slide vertically into the upper part of the cloud and 
allowing the particles to adhere to the slide as a re- 
sult of their own motion; this procedure could conceiv- 
ably have led to a distribution biased toward the smaller 
particles as a result of their greater Brownian motion, 
but we have not made an analysis of that point. Zander 
found the individual particle diameters to be within the 
range ly _< d 3y with the maximum of the distribution at 1.7 
About 1/3 of these single crystals were prisms and most 
of the others were hexagonal plates. Each of these two 
forms showed cavities which resulted in a void space of 
not more than ^ 20%. Zander (1966a) reported that the 
particles which were clustered into aggregates composed 
a minimum fraction of those sampled, but, as will be 
shown below, the smaller number of aggregates were prima- 
rily responsible for the scattering properties of the 
cloud . 

The absolute reflectivity was measured as 78 * 5% at 
X = 1 . 75y and 56 * 5% at X = 2.25y with the numbers being 
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the ratio of the observed intensity to that of a Lambert 
surface with albedo unity. The observed relative reflec- 
tivity from ,95m to 4.0y is shown as the heavy curve in 
Figures 5, 6 and 8. That curve represents the average 
of 5 runs; the individual runs did not depart from the 
average by more than ^5-7% with the largest variations 
being for .95 < X < 1.2 y and 3.0 < X < 4.0y. 

The optical thickness of the cloud in the visible 
region was measured at different times by Zander and 
found to be in the range t = 24 * 3. x(X) is here 
defined such that exp(-i) is the fraction of vertically 
incident radiation which experiences neither absorption 
nor scattering while traversing the cloud. 

The maximum in the observed reflectivity between 
1.5 and 2.0y appears to be at X ^ 1.8y in Zander's 
(1966b) graph but the corresponding minimum in the 
absorbtivity of ice is at X v 1.85y. It is not obvious 
that either the reflectivity measurement or the absor- 
tivity measurement should be in error by as much as AX 
v . 05y, but the discrepancy probably results in part 
from Zander's finite resolving power (AX v . 04y) and the 
fact that the reflectivity falls off more sharply from 
1.85 to 2 . Oy than from 1.85 to 1.5y. In order to have an 
uncluttered comparison of theory and observation we have 
drawn Zander's curve with the relative maximum at X * 1.85y. 



COMPUTATIONAL PROCEDURE 


The Mie scattering computations for the phase func- 
tion (scattering diagram) and the single particle albedo 

w o were made with standard methods for particle size 
distributions specified below. The optical constants 
(i.e., the complex index of refraction, n^ - in^) for 
ice were taken from the tables of Irvine and Pollack 
(1968) which are based on an extensive survey of the lit- 
erature. It is expected that the error in n^ may easily 
be ^ 20%. No temperature correction was applied to the 
tabulated quantities since the tables are for a temper- 
ature (-50°C) sufficiently close to that in Zander's ex- 
periment (~44°C) . The phase function was calculated at 
96 angles [0 ( . 2) 5 (2 . 5) 180° ] . The single particle albedo, 

w , was determined to a high accuracy (±.00001) because 
0 

for albedos near unity the variations in this parameter 

are magnified greatly by the multiple scattering. 

The radiative transfer problem was solved using the 

"double only" method described by Hansen (1968) . At 

each wavelength the phase function was used to determine 

the reflection and transmission functions for a cloud of 

_ 2 5 

optical thickness t v 2 and the corresponding func- 
tions for a layer of twice that thickness were then ob- 
tained from the doubling equations, which are an expres- 
sion of a principle first given by van de Hulst (1963) . 



The process was repeated until an optical thickness 
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t ^ 2 = 128 was reached. 

The phase functions were expanded in cosines of the 
azimuth angle and the first 10 terms in the expansion 
were used; since the incident and emergent angles were 
near normal in Zander's experiment the omitted terms did 
not affect the first three figures of accuracy in the 
results. Hence the primary limitation on the accuracy 
was the number of polar angles used in the Gauss quad- 
ratures; this number was varied and chosen in each case 
to be such that the self consistency (and other checks) 
indicated the errors to be < .1%. 

The computations were made at 25 wavelengths in 
the interval . 95y £ A _< 4y including each wavelength 
in that region at which the absorbtivity has a relative 
maximum or minimum. 

COMPUTATIONAL RESULTS 


The initial computations were made for the particle 

size suggested in Zander's (1966a) original article, i.e., 

the particle diameters were taken to be in the range 

lu < d < 3y with the peak of the distribution at d = 1.7y. 
— — m 

The exact shape of the distribution was unknown but due to 
its narrowness this was unimportant; for the calculations 
it was chosen to follow the "cloud" model distribution of 
Deirmendjian (1964) 
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n(d) a d e-^— Cl) 

m 

where n(d) is the volume concentration at diameter d. 

This size spectrum is shown as distribution 1 in Figure 1, 
the single scattering albedo and the relative cloud optical 
thickness as functions of wavelength are given in Figures 
2 and 3, and some representative phase functions are shown 
in the right half of Figure 4. 

To permit an absolute comparison to Zander's obser- 
vations the theoretical computations were made for 

S(t ; 0 , <j> ; e Q , <t> Q ) 

4 cose cose 

o 

with 0 = 0 = 15° and d> — 4> =0°, 

o of 

where S is the scattering function in the definition of 
Chandrasekhar (1960) . The theoretical reflectivities for 
three optical thicknesses of the cloud are compared to the 
observations in Figure 5. The observed absorption 
features near A = 1.5, 2.0 and 3.0y are present in the 
theoretical curves but the spectral reflectivity does not 
correspond closely to the observations: the major absorp- 
tion features at 1.5 and 2.0y are less than half as 
strong as observed by Zander, the minor features near 
1.0 and 2.5u are entirely absent, and the reflectivity 
is a factor of 3-4 too high in the 3.4 - 4.0u range. The 
agreement would be somewhat improved if the observed curve 
were shifted in the vertical direction but a significant 
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shift is not allowed within the absolute error estimated 
by Zander ( ± 5% ) . 

Perhaps the most likely cause of the discrepancy is 
the size assumed for the particles. Larger particles would 
cause stronger absorption features, and hence, to test this 
effect, calculations were made with particles a factor of 
two larger than those in the first distribution employed. 

The computed reflectivity was changed in the desired sense, 
but to a limited extent, and the fit to the observation was 
still unacceptably poor. 

However , after the above computations were completed 
it was learned from Zander (personal communication) that 
the number of particles clustered into aggregates was 
significant 20% with an aggregate counted as a single 
particle) . This is important since the components of an 
aggregate would not be expected to scatter incoherently 
and hence they should not be treated as individual particles 
in the radiative transfer calculations. The manner in which 
the aggregates should be approximated depends on the ratio 
of the particle size to the wavelength and on the particle 
morphology including the amount of void space. According 
to Zander (personal communication) the distribution of 
aggregates peaked at diameter 9y with some particles being 
as large as 15 y and the aggregates were observed to be 
closely packed with a void space estimated as < 20%. Since 
the particles were tightly packed it may be a reasonable 



approximation to treat them as spheres , even though they 
were not small compared to the wavelength. The diameter of 
the equivalent spheres would be slightly 5%) less than 
the diameter of the aggregates but this fact would probably 
be at least balanced by the fact that Zander found the par- 
ticle diameter to increase with depth into the cloud. There 
fore, calculations were made with a double cloud model size 
distribution (distribution 2, Figure 1) with 20% of the par- 
ticles under the partial distribution peaking at diameter 9 y 
The single scattering albedo and the cloud optical thick- 
ness as a function of wavelength are shown in Figures 2 and 
3. Some representative phase functions are shown in the 
left half of Figure 4 and the cloud reflectivity is given 
in Figure 6. 

The theoretical spectrum for the size distribution in- 
cluding the aggregates agrees quite well with the observa- 
tional reflectivity and the required optical thickness is 
within the range specified by Zander (x = 24 * 3). We note 
the following points: 

The major features in the reflectivity from 1.3 to 
2.2y fit the observations almost exactly. These features 
are very sensitive to the particle size; computations with 
still larger particles [the 'cloud’model with its peak at 
diameter 16y (distribution 3, Figure 1) for which repre- 
sentative phase functions are given in Figure 7j show fea- 
tures (Figure 8) significantly stronger than those for dis- 
tribution 2. The absolute reflectivity with these larger 


particles is too small to be compatible with the observed 


curve . 

The mean extinction diameter is defined by 


r n (r ) Q ext ( r /*) dr 


d = 2 r = 


l 


( 2 ) 


? 

J 

*1 


n(r) Q e xt (r ' A) dr 


where Q t (r) is the extinction cross section for a par- 
ticle of radius r, and rj and r 2 are respectively the lower 
and upper limits of the size distribution. d has the values 
2.16, 10.4 and 20. 8y at X = ,95m for distributions 1, 2 and 
3, respectively. Hence easily measurable variations in the 
reflectivity may occur for changes of less than an order 
of magnitude in the characteristic particle size. 

To further assess the relative contributions of the 
aggregate and individual particles, computations were made 
with a size distribution representing only the aggregates 
(a 'cloud' model with its peak at diameter 9p, di = 2 rj 
=5y, d 2 = 2 r 2 = 1 5 u ) . The results (not shown) were prac- 
tically indistinguishable from those for distribution 2 
(which had 20% aggregates) . The reflectivities in these 
two cases differed by < 3% at each of the 22 common wave- 
lengths at which computations were made, with the reflec- 



tivity slightly lower for the distribution of aggregates 
only. The mean extinction diameter for the latter dis- 
tribution was 11. Op at X = . 95p . These results suggest 
that the reflectivity features depend mainly on the large 
particles in the size distribution as a consequence of 
the fact that those particles largely determine the mean 
particle size. 

The minor absorption feature at 1.03p is hardly vis- 
ible in any of the theoretical spectra and it is unlikely 
that employing nonspherical particles would significantly 
strengthen this depression. Hence it appears that the 
theoretical computations will not yield a feature as strong 
as observed by Zander unless the absorbtivity of ice has a 
greater relative maximum at 1.03y than is indicated by the 
tabulated optical constants. X = lp is in fact near the 
short wavelength limit of reliable determinations of the 
ice absorbtivity because of the high transmission there, 
however, the absorbtivity at 1.03y would have to be approx- 
imately an order of magnitude larger than the value given 
by Irvine and Pollack (1968) in order to account for the 
feature reported by Zander. Perhaps it is more likely that 
the observed reflectivity feature was an instrumental effect 
since, according to W. Plummer (private communication) , the 
wavelength was near the tail of emission of Zander's 
infrared source. Absorption by water vapor probably could 
not have been responsible for the feature because the 
water vapor absorption coefficient is much less at ly than 
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at some longer wavelengths where no absorption was apparent. 

The minor feature at 2.6y is considerably less prom- 
inent in the theoretical spectra than in the observed curve. 
Zander found this feature to vary significantly with tem- 
perature (probably due to changes in the particle size dis- 
tribution) but there is clearly a discrepancy between the 
theory and the observations; this is apparently the result 
of an error in either the tabulated absorbtivity of ice or 
the assumptions employed in our computations. 

In Figure 6 the theoretical reflectivity exceeds that 
observed in the range 3.6 < A < 4.0y, but that is the 
region in which the observations are most uncertain. If 
the deviation is real, though, it could be a consequence 
of the spherical particle approximation since, as mentioned 
by D. G. Rea (personal communication) , aggregates might 
scatter as large particles at long wavelengths and exhibit 
scattering properties characteristic of the smaller com- 
ponents on the aggregate surface at short wavelengths 
where x = - >> 1. Such a scattering behavior has in 

A 

fact been observed for mineral samples in the laboratory 
by Rea, but the possible discrepancy in the cloud reflec- 
tivities must be reconfirmed experimentally before it is 
taken as clear evidence for a breakdown in the spherical 
particle approximation. 

The observed reflectivity of 3-4% from 3.0 to 3.2y 
is perhaps indicative of the noise level in Zander's 
experiment since any particles large enough to yield the 


observed features at 1.5 and 2. Op would have a reflec- 
tivity < 1% for 3.0 < A < 3.2p. 

Our computed results shown in Figures 5, 6 and 8 may be 
compared to the theoretical computations of Sagan and 
Pollack (1967) who employed a two-strea?n approximation for the 
multiple scattering. When their reflectivity curves for 
single particles of given diameters are compared to our curves 
for size distributions with corresponding mean extinction di- 
ameters , it appears that the spectral features are somewhat 
more pronounced in our curves; this, however, is expected 
since their results apply to the intensity integrated over all 
angles of incidence and reflection and not to the intensity at 
the particular angles of incidence and emergence in Zander's 
experiment (which were nearly normal to the cloudtop) . It is 
apparently for this reason that Sagan and Pollack obtain their 
closest agreement to the spectrum of Venus (which is practi- 
cally identical to Zander's laboratory curve for particles 
of diameter ^ 20y while we get the best fit to Zander's lab- 
oratory data with particles of half that size. A preliminary 
comparison of our calculated spherical albedos to their re- 
sults suggests that their approximation is quite accurate for 
the integrated intensity. 


DISCUSSION 


The computed results indicate the sensitivity of the 
near-infrared reflection spectrum to the particle size and 



to the cloud optical thickness; this reemphasizes the po- 
tential use of reflection spectra for cloud identification. 
Our computations were made at rather widely spaced wave- 
lengths and no attempt was made to degrade the theoretical 
spectra to account for the finite resolving power (^.04y) 
in the observations; however, these deficiencies may easily 
be eliminated in cases where more exact computations are 
warranted. 

Due to the uncertainties in the optical constants of 
ice and our incomplete knowledge of the experimental pa- 
rameters such as the particle size distribution, it is 
possible that the close agreement of the major features in 
the theoretical and observational spectra is partly fortu- 
itous. But until further observations and computations are 
made the results must be regarded as evidence that in some 
circumstances randomly oriented nonspherical particles may 
be approximated by spheres; this is a contention which has 
previously been made by Sagan and Pollack (1967) and others. 
However, since the scattering properties of the cloud were 
determined mainly by the aggregates which may have been more 
spherical than the individual plates and prisms, it would 
also be desirable to have reflectivity measurements from 
clouds in which the crystals were not allowed to agglomerate. 

From the minor reflectivity features there is a sug- 
gestion that the optical constants could be in error near 
2 . 6p and perhaps also near ly; hence it would be of interest 
to have the absorbtivities remeasured in those regions. It 



would also be helpful to have the cloud reflection spectra 
reexamined at the same wavelengths. 

Zander's measurements were the first of this type and 
at the time he performed them there were no theoretical cal- 
culations available; it is hence surprising that he obtained 
so many parameters which must be known for an adequate inter- 
pretation (such as the cloud optical thickness, the absolute 
reflectivity, the particle morphology and the size dis- 
tribution) . In future observations it will be important to 
control the experimental conditions as closely as possible, 
especially the particle size distribution which should be 
known at increments of the radius v ly. 

On the whole the agreement between the laboratory and 
theoretical results is sufficiently close to encourage the 
comparison of similar computations to observations of atmos- 
pheric clouds. Measurements from an aircraft of the absolute 
intensity spectra of solar radiation scattered by terrestrial 
ice and water clouds have been reported by Espinola and Blau 
(1965) and Blau et al (1966) in the region 1.15 - 3.6y, 
including a few examples of the angular dependence of the 
scattered light. In their observations the absorption by 
C0 2 and H 2 0 vapor above and within the clouds is significant 
and must be accounted for in a theoretical interpretation; 
we will hence only mention here that the spectral features 
that they observed for cirrus clouds near A =1.5 and 2.0u 
are deep and appear to require particle sizes^at least lOy 



in diameter. Further observations including the angular 
pattern of the scattering, especially in the region 1 < X 
< 2.5u where gaseous absorption is small, would be helpful 

for making a complete analysis. 

Finally, the results of this paper also suggest that 
it is possible to obtain meaningful conclusions from the 
theoretical interpretations of the near infrared reflectivity 
of Venus (Sagan and Pollack, 1967; Hansen and Cheyney, 1968). 
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FIGURES 


Fig. 1. Cloud particle size distributions normalized to 
unity at their maximum points. The single and multiple 
scattering results in Figs. 2-8 were made for these three 
size distributions of spherical ice particles. Distribu- 
tions 1 and 2 coincide for diameters less than 3m . 

Fig. 2. Single scattering albedo (u> ) for the three 
particle size distributions shown in Fig. 1. The wave- 
lengths at which the single and multiple scattering compu- 
tations were made are indicated by dots on the curve for 
distribution 3. 

Fig. 3. Cloud optical thickness normalized to unity at 
A = ,95m for the three particle size distributions shown 
in Fig. 1. The interpolated curves between the wavelengths 
of the computations (indicated on curve 3) may not be accurate 
for 2.8m < A < 3.2m, but that is unimportant for the mul- 
tiple scattering results (Figs. 5, 6, 8) since the inten- 
sity is so low in that wavelength region. 

Fig. 4. Single scattering phase functions (scattering 
diagrams) normalized to unity, at representative wave- 
lengths in the near-infrared for two of the particle size 
distributions shown in Fig. 1. The vertical scale for 
each successive wavelength is displaced by an order of 
magnitude; the wavelengths to which the numbers on the 


vertical scale apply are indicated in parenthesis. 

Fig. 5. Near-infrared reflectivity of an ice cloud. The 
solid curve was observed in the laboratory with a resolu- 
tion ^ . 04y by Zander. The theoretical curves are for the 
particle size spectrum without aggregates (distribution 1, 
mean extinction diameter v 2.16u). The optical thickness 
applies to A = . 9 5 p . The reflectivity is defined as the 
ratio of the observed or theoretical intensity at the 
angles in Zander's experiment to the intensity from a 
perfect Lambert surface at the same angles. 

Fig. 6. The same as Fig. 5 with the theoretical curves 
for the particle size spectrum including aggregates 
(distribution 2, mean extinction diameter 'v 10. 4y). 

Fig. 7. Same as Fig. 4 for size distribution 3. 

Fig. 8. The same as Fig. 5 with the theoretical curves 
for the particle size spectrum of very large particles 
(distribution 3, mean extinction diameter ^ 20. 8y). For 
A > 2 . 8y the three theoretical curves are indistinguishable. 
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